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We present a simple argument to account for crystallization of hard spheres under the action of 
a gravitational field. The paper attempts to bridge the gap between two communities of scientists, 
one working on granular materials and the other on inhomogeneous liquid state theory. 
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I. INTRODUCTION 



Sedimentation of colloidal particles in a gravitational field has been the subject of continued interest since the 
early work of Jean Perrin. The external gravitational field breaks the translational symmetry, producing an altitude 
dependent density profile. For a dilute suspension in which the interactions between the colloidal particles can be 
neglected, the density profile takes a simple exponentially decreasing form. In fact it was this observation which has 
first allowed Perrin in 1910 to obtain an estimate of the Avogadro number and the Boltzmann constantly. This early 
' work has stimulated much of the modern development of inhomogeneous fluid theory. 

A few years ago it has been observed in Monte Carlo simulations that a system of hard spheres can become unstable 
under the action of a sufficiently strong gravitational fieldcl. This instability appeared as a sudden crystallization of the 
bottom layers of the sample confined to a tree dimensional rectangular box. In order to understand this transition the 
authors appealed to the Weighted Density Functional theory (WDF). Unfortunately the WDF is quite computationally 
demanding and is not very physically transparent. Perhaps it is for this reason that this work went unnoticed by the 
community of physicists studying granular materialsti Thus, the fluid-solid condensation under the gravitational field 
has been recently rediscovered in the context of granular materials. Here the instability appeared as the break down 



of the Enskog kinetic theory for elastic hard spheresa, which Hong associated with the formation of a solid layer at 
the bottom of a sampleQ. We must stress, however, that the identification of the break down point of the Enskog 
equation with crystallization is far from obvious. In particular, .it is possible to show that at equilibrium the Enskog 
■ theory is equivalent to the Local Density Approximation (LDA)cl. Within the liquid state theory it is well known that 
the LDA cannot be used to obtain reliable information about the fluid-solid coexistence. The reasosufor this is that 
the LDA fails to account for strong density variations in a highly structured phase, such as the solidOu. Furthermore, 
in the absence of an external field, the LDA for hard spheres reduces to the usual Carnahang-Starling equation of 
stated, which carries no information about the entropy driven crystallization of a homogeneous hard sphere fluid. It is 
difficult, therefore, to understand how a theory which knows absolutely nothing about the existence of an instability in 
a homogeneous system, suddenly develops this "knowledge" in a presence of a gravitational field. Thus, it is dubious 
that the break down of the LDA — and, therefore, of the Enskog theory — can be associated with a phase transition, 
q | and not with the failure of the approximations underlying the theory. 

In view of all this, it seems worth while to present a simple argument to account for crystallization of hard spheres 
under the action of an external gravitational field. 
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H II. LATTICE GAS IN A GRAVITATIONAL FIELD 

Before entering into the discussion of a .-phase transition in the system of hard spheres, we shall explore a much 
simpler model of fluid and granular mattem in a gravitational field — a non-interacting lattice gas. Consider a three 
dimensional rectangular volume V — L 2 H, divided into cubic cells of volume a 3 . Each cell can be occupied by at 
most one "cubic" particle of mass m and volume a 3 equal to that of a unit cell. The constant gravitational field of 
strength g acts in the — z direction. Since the particles do not interact, in the thermodynamic limit, it is a simple 
matter to write an exact Hclmholtz free energy functional, 

^^^=^{[\-a{z)]\n[l-a{z)]+a{z)\na{z) + Pmagza{z)} , (1) 
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where /? = l/fc^T, and u{z) is the average occupation of the row z. The equilibrium density profile can be obtained 
by minimizing the functional (]l|), subject to the constraint of conservation of the total number of particles N, 

oo 

£V(z) =Na 2 /L 2 = n. (2) 

Performing a simple calculation we find, 

where the chemical potential [i is the Lagrange multiplier associated with the constraint (|^) . It is convenient to define 
the reduced parameters, f3* = (3mga and ji* = u/mga. In the limit —0fi 3> 1 the reduced chemical potential is very 
well approximated by u* = —n. The density profile now takes a particularly simple form, 

ml 

' (z) = ^l + e0-(»-»> • (4) 
On the other hand the static pressure inside the lattice gas satisfies, 

d_P{z} 

—, = -9ap(z) , (5) 
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where to simplify the calculations we have passed to the continuum limit. Eq. ^| can now be integrated yielding, 
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We note that the pressure is a uniformly decreasing function of the height z. At high altitudes — where the average 
density is small and the excluded volume constraint is irrelevant — the density profile reduces to the exponentially 
decreasing form first obtained by Perrin. 

III. THE HARD SPHERE FLUID IN THE GRAVITATIONAL FIELD 

It is evident from the analytic expressions derived above that the lattice gas — in spite of its hard core repulsion 
expressed through the constraint of one particle per lattice site — is completely stable for all temperatures. It is 
reasonable to suppose that to zeroth order the density profile of a hard sphere fluid in a gravitational field should 
also be well approximated by the Eq. (Q). In fact it has been found that the Fermi- like distribution (Q), provides an 
excellent fit for various granular materials including the hard spheres&Ej. 

Unlike the non-interactpg lattice gas, the homogeneous system of hard spheres is known to undergo a struc- 
tural fluid-solid transitionU when the bulk pressure reaches Pfr— p 14fceT/a 3 . This entropic transition was first 
observed in Moste-pCarlo simulations almost half a century agcO, and has since been confirmed by various theo- 
retical methodgl3'H , lij. In the presence of a gravitational field the pressure inside the suspension varies uniformly 
with maximum value located at z = 0. In the limit (3*n 3> 1, the pressure at the bottom is well approximated by 
P(0) = mgn/a 2 . Therefore, the first layer will crystallize when, 

= ragn ^ Uk^T ^ 
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or when 
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This is a discontinuous first order transition. If the temperature is lowered further, additional fluid layers will solidify 
and a growing crystal will coexist with a diminishing fluid phase. The number of layers in the crystal can be determined 
by comparing P(z) with Pi,. We note that in the region of phase transition, 

(3*n w -0V ~ 14 > 1 , (9) 
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so that all the approximations adopted above are indeed justified. 

We have presented a simple argument for crystallization of hard spheres in a gravitational field. This transition is 
nothing more than a realization of the well known entropically driven liquid-solid phase separation of a homogeneous 
hard sphere fluid. The presence of a gravitational field leads to an altitude dependent pressure, which increases with 
the depth of suspension. If the temperature is reduced bellow the critical value T c , the pressure on the bottom layers 
becomes so large as to lead to their crystallization. 

A comment should be made about the real granular materials. For macroscopic particles, temperature is an 
irrelevant parameter. Thus, in order to study dynamic equilibrium or steady-state, kinetic energy must be supplied to 
the system through, say, a vibrating bedE2l. If the collisions between the particles are elastic, in the disordered (fluid) 
state we can think of the kinetic energy as a form of "temperature" . Appealing to the same argument as above, we 
reach the conclusion that the granular materials also undergo phase separation under the action of a gravitational 
field. Some care, however, must be taken with this argument. It is well known that unlike for fluids, the pressure 
inside a container filled with packed granular matter saturates at some valueEJ. The excess pressure is transfered to 
the walls and is equilibrated by friction forces. At high packing fractions, where the solidification takes place, the 
stress is not uniform but is concentrated along force lines. This unusual characteristics can significantly affect the 
location of the phase transition in real granular materials. 
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